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In the one-loop approximation we derive the equation of motion for a classical scalar field 
<f c (t) with the back reaction of particle production included. Renormalization of mass and 
qq | couplings of ip c is done explicitly. The equation is non-local in time, but can easily be treated 

■ perturbatively or numerically. For the weak trilinear coupling of the external field to the pro- 

duced particles, the new equation gives the same solution as the familiar one with the Tip c 
term. For a stronger coupling and other types of couplings the results are significantly different. 
The equation can be applied to the universe heating by the inflaton decay and to spontaneous 
baryogenesis. 



X 



PACS: 98.80.Cq, 14.80.-j, 95.30.Cq 

1 Introduction 



In the classical theory of particle production (see e.g. [|J) the external time-dependent field is 
assumed to be "fixed", i.e. not affected by the back reaction of the produced particles. Though it 
is a very good approximation for the situation that can be realized under laboratory conditions, it 
is not so in cosmology, when an external cosmic field decays and transfers energy to the produced 
particles. Usually the back reaction is described by an addition of the term Ttp c into the usual Klein- 
Gordon equation governing the evolution of a (spatially homogeneous, though not necessarily so) 
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classical scalar field: 

<p e + 3Hp c + V'(<p e ) = -Tip, , (1) 

where H = a /a is the Hubble parameter describing the expansion of the universe (this term is 
absent in flat space-time), V((p c ) is the potential of the field ip c , the prime denotes the derivative 
with respect to the field cp c , and T > is the decay width. One argument in favour of a modeling 
the back reaction in this way is that the solution of eq. (|l|), in the case of the simplest harmonic 
potential, V((p c ) = m 2 ip^/2, and for the trilinear coupling to the produced particles, Ci n t ~ <£cX*Xi 
has the form: 

(p e (t) = exp(-n/2) cos (y m 2 -r 2 /4 tj , (2) 

which correctly describes the behaviour of (p c known from perturbation theory (it is assumed here 
and in what follows that the space-time is flat). However, it is easy to see that already in the 
case of more complicated potentials and couplings, e.g. for V(<p c ) = X^p^/A and/or Ci n t ~ <PcX*X, 
this ansatz is not applicable. Moreover, different equations may have the same solutions and the 
coincidence of the solutions in special cases does not imply the equivalence of the theories, as we 
will see below in some examples. 

In what follows we will derive the equation of motion for the classical homogeneous scalar field, 
(p c (t), for different kinds of interactions with the produced quantum particles and different forms 
of the potential, V(<p). In fact, the back reaction term does not depend on the form of V(cp), but 
is determined only by the form of the interaction of (p c with the produced particles. On the other 
hand, the particle production rate is sensitive to the form of the potential V(ip). The calculations 
here are made in flat space-time, but it is straightforward to generalize them [[| to the case of 
the Friedman-Robertson- Walker metric, because the latter is conformally flat. The results can be 
applied to the universe heating by the decay of the inflaton field and to spontaneous baryogenesis. 
Our approach here is very close to that of ref . 0] . 
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2 Trilinear scalar coupling 



Let us first consider the simplest case of the trilinear coupling of a scalar field ip (which below will 
be taken as a classical spatially homogeneous field, (p c (t)) to the massless scalar quantum field x- 
The Lagrangian of the fields (p and x nas t ne form: 

C(<p, x) = \{d^) 2 - V{y) + \{dx? + f^PX 2 ■ (3) 
The corresponding exact quantum operator equations of motion can be written as: 

<p-A<p + V\<p) =f X \ (4) 
d 2 X = 2f<p X - (5) 

We will assume that there exists a classical homogeneous field <p c (t) = (<p), where the brackets 
mean averaging over quantum vacuum state of operators ip and x- The equation of motion for ip c 
can be obtained by the quantum averaging of eq. (|j) in the presence of the classical field tp c [§]: 

^ + V\ V )) = (f X 2 ) ■ (6) 

Taking average of the l.h.s. is trivial, it reduces to the substitution {ip) = ip c (t), and averaging of 
the r.h.s. can be done using the formal solution of eq. (|5|): 

X (x) = xo(x) + 2/ J d 4 y G R (x, y) ^ c {y)x(y) , (7) 

where G R {x,y) is the retarded Green's function and X o is the free field operator. The latter is 
canonically quantized according to: 



(2vr) 3 / 2 v^ 

where a? and at are annihilation and creation operators of momentum k, obeying the commutation 



relations 
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(2^) 3 2A;o 5 3 (k - k'). 



The contribution of expression (^) into the r.h.s. of eq. (|6|) gives zero result to the first order in 
/. To be more precise, the condition of vanishing of (xx) to the first order in / corresponds to 



elimination of tadpole diagrams and is achieved by a shift of ip c . To the second order in / the 
contribution is given by: 

HXX) = 4/ 2 / I d A ye- ik ^G R (x,y)My) • (9) 



(2ir) 3 2k 

Using the retarded Greens function for massless scalars: 

Gf(r) = ^- r S(r-t), (10) 

we can make almost all the integrations in eq. @ and obtain the equation of motion of the classical 
scalar field, (p c , with the r.h.s. which describes the back reaction of the produced quanta of scalar 
field, X-, on the evolution of the field ip c : 

<Pc + V'(<p c ) = ^J Q -^(t-r), (11) 

where t{ n is an initial time, when the particle production was switched on (it is assumed that 
t > ti n ). The result is valid for the coupling of <p c to x given by expression @. For other forms 
of couplings and for produced particles with non-zero spin the r.h.s. would have different forms 
(see below). The term describing the back reaction of the produced particles is determined by the 
form of the interaction Lagrangian of the field (p c with the produced quantum fields, and does not 
depend upon the potential V(ip c ). This is in contrast to the naive description of particle production 
by the term Tip c with a constant T, suggested in earlier literature. It may approximately describe 
the realistic situation only for the harmonic potential, V = m 2 (p 2 /2. If one wants to mimic particle 
production for other types of the potential by a similar local term, one would immediately conclude 
that such a term must have different forms for different potentials. 



The integral in the r.h.s. of eq. (11) is logarithmically divergent at r — > 0. This divergence is 
related to the mass renormalization and can be regularized in the following way. Let us separate 
the integral in the r.h.s. into two parts, one from to some t\ and the other from t\ to (t — ti n ). 
In the first integral let us subtract and add <p c (t). We obtain: 

<p c + V'(cp c ) - In — 

= T=2 T^-r)-^)] + ^ ?Pe(t-r). (12) 



4tt 2 Jo t An 2 Jt x r 



The integrals in the r.h.s. are now finite. The logarithmically infinite contribution in the l.h.s., 
related to the ultraviolet cut-off e —* 0, is taken out by the mass renormalization, so with a possible 
bare mass term in the potential, V mo = m\ip 2 c j2 (here m is the bare mass of ip c ), we obtain: 

ypc + m 2 (h)^ c + [V'(<p c ) - V^(tpc)] 

■i-S _ - [<p e (t - t) - <p e (t)] + ±3 / -Vc(t-r), (13) 



47T 2 Jo T 47T 2 r 

where t\ is an arbitrary normalization point and the "running" mass is m 2 (t\) = mPfe) - 
(/ 2 /47r 2 ) h\{t\/t2). In eq. ( |i~3| ) we explicitly separated the massive part, V mo , in the potential 
so that the term in the square brackets vanishes for the harmonic potential, V(ip) = mg(^ 2 /2 



In the limit of a small coupling, /, and for the harmonic potential, equation ( p~3[ ) can be solved 
analytically 0. In this limit its solution coincides with that of the naive eq. (0). We will look for 
the solution in the form: 

ip c = A(t) • sin(mi) , (14) 



where A is a slowly varying function of t. Substituting this into eq. fllq) we obtain: 

— 2Amcos(mt) = ^— ^ — Acos(mt) , (15) 

and correspondingly A(t) ~ exp(— Tt/2) with T = / 2 /87rm. It is the correct decay width of 
ip c — > 2\- Thus we found, that in the limit of a weak coupling both equations (|l|) and (|i~3| ) have 
the same solution, though the equations themselves are quite different. 

For the case of the harmonic potential one can reduce equation (|ll|) to a more familiar form by 
making the Fourier transform: 

— exp(-iu;t)tp c (uj) , (16) 
-oo ^vr 

where the spectral function (p c {u) should be analytic in the complex upper half-plane of oj to ensure 
vanishing of <p c for t < ti n (below we take tj n = 0); if <p c = const ^ for t < 0, this constant value 
can simply be subtracted from <p c . Assuming that t\ is sufficiently small {rati <C 1), so that the 
integral from to t\ in the r.h.s. of eq. ( p^ ) or (|l3|) can be neglected, we obtain: 

/ f2 + \ f2 r oo pt A T 

ipM [ral-±-An±-u? = / cfte"-*/ - ^ c (t - r) . (17) 



47T 2 £ / 47T 2 J-oo Jh T 



For a large t the integral over dr depends weakly upon the upper limit, so that it can be pushed 
to infinity and the integration over t can be done explicitly: 

p2 ± \ f2 



(p c (u) ml 



K m U = Ktfic / — [COS LOT + ism LOT) . (18) 

A-k 1 e I 4tt z Ju t 



The second integral in the r.h.s. is equal to i mT sign(cc>)<y5 c (u;), while the first one, logarithmically 
divergent in the lower limit, gives mass renormalization in the momentum space. With the account 
of this term the renormalized mass can be written as: 

f2 h _ f_ r°° dy * 2 

e 47r 2 Jutx y 
where A = 1/e is the ultraviolet cut-off. Finally we obtain: 



2 , *i f 2 [°°dy 2 f* a 2 

m n — — ^ In —^r / — cos y = m n — — ~ In h (nnite terms) = m^.^ , (19) 

u 4-7T 2 e 47r 2 Jut, w u 4-7T 2 a; 



<^ c (a;) m 2 en - Lo 2 - i m ren T sign(w) =0. (20) 

This is a simple linear equation in momentum space. A similar equation for an unstable particle 
moving in an external potential (gravitational field) has been derived in ref. || . Eq. (^) is different 
from eq. (||) with the T^c-term, but in the limit of small T/m it has the same solution (H). With an 
increasing T/m the solution of eq. (||) ceases to oscillate and turns into an exponentially decaying 
one. For an even larger ratio T/m, especially in the unphysical limit r>m, the solution of eq. 
(|l]) tends to a constant. The solutions of eq. (^) behaves differently. For any value of the ratio 
T/m the solution of this equation has both oscillating and exponentially decreasing factors. This 
behavior is demonstrated by the numerical solution of eq. (|l3l), see figs. 1. Recall that eq. (|l3|) is 
equivalent to the algebraic one (^0|) for the Fourier transformed field amplitude, only in the limit 
of a large mt. In the general case eq. (12) or (13) should be solved numerically. 

There is a subtle point related to the definition of the value of the renormalized mass. It is clear 
that m ren in eq. ( p0|) and m in eq. @ (correspondingly m in V(ip c ) in eq. (|l|)) are normalized at 
the same point and should be taken equal when one compares the solution of these equations. On 
the other hand, the value m 2 (t\) in eq. (|i~3|) is different from m 2 en , because of a finite contribution 
to m 2 from the r.h.s. of this equation, and depends therefore on /. This explains the presence of 



the coefficient a ^ 1 in the envelope exp(— aTt/2) to the numerical solution of eq. (13), presented 
in figs. 1. 
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There is no simple way to find an analytic solution for more complicated potentials, when the 
non-perturbed equation for <p c (t) is non-linear. For the case when the potential is dominated by 
the quartic term, V(<p) = X^ip^/A, the equation can be written as: 

z » + z * = h 2 f ^z(x - y) . (21) 

Jx! y 

Here the effective coupling constant is h? = f 2 / (47r 2 A^q) ■ We introduced the dimensionless 
quantities x = \/A~4^o* and z(x) = ip/ifo, where ipo is the characteristic magnitude of the field 
if. Prime means differentiation with respect to x. The contribution to the integral from to x\ 
is neglected as in the previous case (but it can easily be taken into account). We will solve this 
equation for ip(t = 0) = ip and ip{t = 0) = 0. 

In the limit of h = this equation is solved in terms of Jacobi elliptic functions (see e.g. j|]): 

z {x) = cn(x,l/v / 2) • (22) 



In what follows we will use the notation cn(x, 1/V2) = cn(x). It is easy to find the first order 
corrections in h 2 in the limit of relatively small time, h 2 x <C 1. We will look for the solution in the 
form: 

z(x) = A(x) cn \J3(x)] , (23) 

where A{x) is a slowly varying function of x and (5'{x) = A(x), which ensures vanishing of the 
lowest order terms in the l.h.s. of eq. (21). The function A(x) differs from the one used above 
by the constant factor ip^. Substituting expression ( |23| ) into eq. ( |2T| ) and neglecting the term 
proportional to A" we obtain: 

3A'A-^ [cn (5{x)} = h 2 f ^-A(x - y)cn [(3(x - y)} . (24) 

In the first order in h 2 x the solution is straightforward: 

A(x) « 1 - 0.62/i 2 x ~ 1 - 0.016/ 2 f \/%<p ) 1 1 . (25) 



Unfortunately we failed to find an analytic approximation to the asymptotic decrease of the am- 
plitude of the oscillating field if at very large x, which characterizes the particle production rate at 
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a large time. The problem is that A{x) — * for x — > oo and (3{x) — > const. Correspondingly the 
decomposition of z as a product of quickly and slowly varying factors becomes invalid. However, 
the numerical treatment of the problem is quite simple and straightforward. The r.h.s. of equation 



(21) gives a nonzero contribution into the mass of the field 92, as is easy to see a negative one, 
5m? < 0. Correspondingly the solution, found numerically, asymptotically tends to the nonzero 
value, v'jin = —5m 2 /X. A counter-term 5m 2 ip 2 /2 can be added to the potential, V(<p), to ensure 
vanishing of the renormalized mass. The equation is solved in this case too and, as expected, ip — * 0. 
The decrease of its amplitude at small h?x agrees quite well with the approximate result ( p5[ ) as is 
seen in figs. 2. 

In connection with the derivation of the equation of motion for the primary field <p, given above 
and in the following sections, a very important question may arise. We derived equation ( |l3| ) 
perturbatively including only one-loop contribution into {xx)- Could one go beyond perturbation 
theory using this equation? We do not have a rigorous answer to this question (see also the 
discussion in the last section). However, there are quite many examples, that an equation with a 
perturbative potential (for example the Schroedinger equation with Coulomb potential) permits to 
go beyond perturbation theory. Moreover, the fact that equations (|l|) and (|l^) are different even 
in perturbation theory have important implications ||, |j| e.g. for spontaneous baryogenesis || 0. 



3 Quart ic scalar coupling 

The calculations are essentially the same for other types of coupling of ip c to quantum fields. For 
the interaction of the form Ci n t = \2f\x 2 the equation of motion for ip c takes the form: 

<p c + V '(<p c ) = £L Vc{t ) J q — £(t - r) . (26) 

The logarithmic divergence at r — > can be removed by renormalization of the coupling constant, 
A4, in the self-interaction potential, V\(tp c ) = A4^/4. It can be done exactly in the same way as 
it was done above for mass renormalization. We separate the integral in the r.h.s. of eq. (^) into 
two parts: from to t\ and from t\ to (t — tj n ). In the first integral we subtract and add <p 2 (t) and 
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obtain: 




<p 2 c (t-T)-<p 2 c (t) 



— V 2 c(t-T) , 



(27) 



where \4(ti) is the coupling constant of the ip 4 self-interaction, renormalized at the normalization 
point ti, by the loop with the field %, A4(ii) = A4(i2) — (Al/^vr 2 ) ln(ti/t2). The integrals in the 
r.h.s. are now finite. We assume that t\ is small, mt\ < 1, while the solutions will be taken in 
the limit mt > 1. Here we neglect the loop with the field ip itself or, in other words, neglect the 
self-production of 92-quanta by the classical field (p c . This is physically sensible if A4 < A2. 

A perturbative solution in this case is slightly more complicated than that given in the previous 
section. We assume that the potential is dominated by the harmonic part, V(ip c ) « m 2 ip 2 ./2, and 
that the coupling constants are small, A4 < A2 < 1. To the lowest order in the coupling constants 
the solution has the form: 

ip c (t) = A(t) cos(mt + $) + A 3 (t) [at sm(mt + + c 2 sin(3mt + 3$) + c 3 cos(3mt + 3$)] , (28) 

where $ is a constant phase determined by initial conditions. The slowly varying amplitude A(t) 
and the coefficients Cj are given by: 




I-1/2 



A(t) =A in 1 + 



167T m 



(t tin) 



(29) 




(30) 



c 2 = -A^/(128vrm 2 ) , 



(31) 



c 3 = A/(32m 2 ) , 



(32) 



where: 
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The term proportional to c\ in equation ( |2S| ) rises linearly with time (in addition c\ itself rises 
as In mt). This linear rise is the well known resonance behavior associated with the equality of 
the eigen-frequency of the non-perturbed equation and the perturbative force, coming both from 
the term proportional to A^j? and from the terms in the r.h.s. of eq. (27). Moreover, if A4 is 
non-negligible (to be more precise, the combination X^A^/m 2 ), then the term X&ipqipi (where (fo 
is the zero-order approximation and ipi is the first-order correction) would induce a parametric 
resonance and an exponential rise of 99. Of course these resonances do not have physical sense, 
because the original equation does not possess any instability. The "resonances" (both the usual 
and the parametric ones) are quickly turned off with rising time and the large time behavior remains 
non-resonant, as is confirmed by the numerical calculations. A resonance amplification of the boson 
production rate would take place if the accumulation of the produced bosons is taken into account. 
To do that one has to make quantum averaging of the equations of motion, not over the vacuum 
state, as is done here, but over the states with non-zero boson occupation number. It will be done 
in a subsequent work Q. 

The amplitude of the field, ip c {t), decreases in accordance with expression (p9|). This decrease is 
much slower than the exponential decay of (p c (t) considered in the previous section. Correspondingly 
the efficiency of particle production in the case of the quartic coupling is much weaker than that 
induced by the trilinear coupling in the previous section. The usual description of the back reaction 
of particle production by the Ytp c term in this case is very far from reality. Correspondingly the 
approach of the earlier papers |9|, [l(J (see also the recent review [11]) should be modified. The 



results of the numerical solution of eq. ( p7| ) are presented in figs. 3. For a weak coupling, i.e. for 
X\w? /2-7T 2 Af n <C 1, they agree quite well with the perturbative result (|28|). The fitting coefficient 
a l appears because of renormalization of the parameters entering the l.h.s. of eq. ( |27| ) due to 
the interaction term in the r.h.s., while the dependence on time remains the same as that given by 
eq. (H). 

The case of the quartic potential, V((p) = X^^/A, can be treated similarly to that in the previous 
section. The solution is searched for in the form (p3), so that the l.h.s. of eq.(pl) remains the same, 
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but the r.h.s. becomes different: 



SA'A^cn [{3(x)\ = -^-cn [(3(x)} f ^A 2 (x - y)cn 2 (J3(x - y)) , (33) 



dp n 2vr 2 A 4 J y 

where x = y/X^ip t. 



The decrease of the amplitude of the oscillations is described by the expression: 



1 

1 + 0.0085(A!/A4)c / ?o£ 



A (t) « t , nnnoea2/ , — • (34) 



To obtain this result, only the terms cos(px) and cos(3px) (with p ~ 0.85) in the Fourier decom- 
position of cn[/3(x)] (see e.g. Q) were taken into account. Again this behavior agrees quite well 
with the numerical solution for small couplings, see figs. 4. As above the fitting parameter a / 1 
is induced by renormalization. 

It is interesting to note, that for the quartic coupling considered in this section the oscillations 
asymptotically decrease faster in the case of the (/^-potential, eq. (p4|), than in the case of the 
</? 2 -potential, eq. (]29|). Further, comparing figs. 1 and 2 with figs. 3 and 4 (and the approximate 
behaviour of the envelopes) one realizes, that for the trilinear coupling, f^px 2 , the scalar field, <p, 
disappears much faster than for the quartic coupling. This is what one would expect, since the 
probability of the reaction induced by the coupling \2 i f 2 X 2 is quadratic in the <p- number density, 
whereas for the trilinear coupling the probability is linear in the ^-number density. 

4 Coupling to fermions 

When the classical scalar field couples to quantum fermions the calculations are slightly more 
complicated, but still straightforward. We will consider the Lagrangian density where the classical 
scalar field, ip c , couples to the spin 1/2 field i/j through the Yukawa coupling: 

where tp is quantized according to: 

= J { jll/2E £ ( a > S ^ e ~ iPX + b P vS (Py PX ) . ( 35 ) 
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where a and are annihilation and creation operators at momentum p and spin s, obeying the 
anti-commutation relations: 

= {b r p ,bf} = (2vr) 3 5 3 (p- q)S rs . 
Summation over spins, which we need in what follows, is achieved with the usual relations: 

^ u s (p)u s (p) = "ft + m and v s (p)v s (p) = f — m. 
The equations of motion are: 

<p-A<p + VX<p)=-gW>, (36) 
ip^) + gipip = , (37) 

-ta M (^y)+^ = o. (38) 



In the background of the classical field (f c (t) eq. (p7|) is formally solved by: 

VKx) = i> {x) +ig J d A y Gf(x, y) i>(y)<p c {y) , (39) 

where ipo is an initially free fermion field and the retarded Green's function for fermions is given 
by: 

Here we will consider the case of massless fermions, which is technically simpler, and take fermionic 
masses into account in a subsequent paper Q. There are some subtleties in the calculations in 
comparison with the case considered above of production of massless scalars, especially regarding 
renormalization and regularization. However, the calculations can essentially be reduced to the 
scalar case. The fermionic Green's function is expressed through the scalar one as: 

Gf(x,y)=ip x G?(x-y) . (40) 



The vacuum expectation value of the field operators in the r.h.s. of eq. ( p6| ) is equal to: 

9${x)4{x)) =g 2 J d 4 yip c (y) dG ^ X x ~ y) Tr(M^h,Mx)) ■ (41) 

12 



Using the representation (35) and taking vacuum expectation value and spinor trace we obtain: 

dG?(x - y) dS(x - y) 



g(ip(x)Tp(x)) = Ag 2 / d 4 yip c (y) 



where: 



S(x - y) 



d 3 k 



dx^ 



a ik(x-y) 



dx^ 



(42) 



(43) 



(27r) 3 2£ 

is the same function (coming from vacuum averaging of free field operators) which enters the r.h.s. 



in the scalar case. It satisfies the condition: 



d 2 S = 



(44) 



while the Green's function satisfies: 



d 2 G{x-y) = 5 A {x-y) 



(45) 



Using these two equations and the identity: 



2duSduG = d 2 (SG) - Gd 2 S - Sd 2 G 



(46) 



we obtain: 



<pc + V'(ip c 
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2 d 2 /"*-*»»• dr 



ip c (t - t) - g 2 A 2 tp c . 



(47) 



4vr 2 dt 2 jo 

The last term, proportional to the ultraviolet cut-off parameter, A, renormalizes the mass of ip c . The 
term under the integral sign is the same as the one in the scalar case eq. (|Tl|) and its logarithmic 
divergence at r = can be treated in a similar way. However, instead of the term </? c ln(ii/e) 
we obtain (p c ln(ti/e). This corresponds to wave function renormalization (renormalization of the 
kinetic term) at the time moment t\. Finally we obtain: 



9 2 d 2 



+ V'. en {ip c ) 



*! dr , . , , N1 /"* tin dr , , 
— [<Pc(t - T ) - <Pc(t)\ + / — <Pc(t - < 



4n 2 dt 2 {Jo t 

where V ren (ip c ) is the potential of ip c with renormalized mass. 



(48) 
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The perturbative solution in the case of a harmonic potential, V(tp c ) = m 2 ip 2 /2, is essentially the 
same as in section ^. If we assume (p c = A(t)sm(mt), where A(t) is a slowly varying function of 
t, we find that A/A = g 2 m/(167r). For a small T it is essentially the same result as can be found 
from the naive ansatz with the Tip c term in the r.h.s. of eq. (|l]) using ^ = g 2 m LPc / Q. 
For the harmonic potential, V(ip) = m 2 ip 2 /2, and the Yukakwa coupling to fermions, giptpip, the 
particle production is quite similar to that considered in the first part of section 2. The results of 
the numerical calculations for this case are presented in figs. 5. 

For the quartic potential, V(ip) = A4</9 4 /4, we will look for the perturbative solution in the form 
(|23|), In direct analogy with the previous sections it is easy to find, that the amplitude A(x) decays 
as: 

A(t) « 1 ^ nm f° 2 ^ - ■ (49) 
The results of the numerical calculations are presented in figs. 6. 

Now, one could also consider non-renormalizable couplings, such as Ci n t = Kip 2 ipip, where k has 
dimension (mass) -1 . The corresponding equation of motion has the form: 



Vc + V'fcftO = ~ <p(t) ^ j Q "" ^ <p e (t - t) 2 . (50) 

An approximate analytic solution for the potential V(<p) = m 2 ip 2 /2 can be found along the same 
lines as presented in the beginning of section 3. The essential part of the solution behaves as: 

( V 2 K 2 \ " 1/2 

ip(t) = cp I 1 + — mtj cosmt . (51) 
This result describes the numerical solution quite well (see figs. 7). 

If the potential is dominated by the A4<£> 4 -term, eq. ( p0[ ) can be solved with the ansatz (p3|). The 
amplitude decays as: 

1 /3 

A(t) = (l + 0.075y%K 2 pit) ~ (52) 



2 



d 2 /"* tin dr 



The decay is astonishingly slow. It can be compared with the numerical solutions presented in 
figs. 8. The agreement is quite good if the above mentioned effect of renormalization is taken into 
account. 
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The numerical solution of eq. ( |50|) can be compared to the solution of eq. (48) and one sees that, 
as in the bosonic case, the scalar field, <p(t), falls off much faster for the coupling gipipip than for 
Ktp^'ipif}. The reason being, that the latter goes quadratically with the number density. 

5 Discussion and conclusion 

We have seen that with a weakly coupled external field, <p c (t), the effect of particle production 
on its evolution can be described by the addition of an extra term into the usual Klein-Gordon 
equation. This extra term is given by an integral over all previous history of the field, ip c (t), arid 
the concrete form of the integrand is determined by the form of the interaction of ip c with the 
produced particles (see examples considered above). In all the cases it has been assumed that the 
interaction is bilinear in terms of the produced quantum fields: 

£ int = F b ( i p)x*X + F f (<p)^ , (53) 

where x an d ip are bosonic and fermionic quantum fields respectively. The calculations have been 
done in the limit of a negligible mass of the produced particles and in flat space-time. Lifting both 
restrictions is straightforward but tedious and will be considered elsewhere ||. 

The most essential restriction is the weakness of the field, tp c (t), such that the Green's functions 
of the produced particles can be taken as the free field ones. This assumption is justified if the 
effective mass of the produced quanta, generated by the interaction (53), is small in comparison 



with their characteristic energy. The latter is essentially the frequency of the oscillations of <p c (t). 
In this approximation our equation adequately describes the decrease of the amplitude of <p c {t) 
induced by the particle production. 

Another simplifying assumption, which was made above, is that the number density of the pro- 
duced particles is small, so that neither Fermi-Dirac suppression nor Bose-Einstein enhancement is 
essential. Correspondingly the production of fermions would be further suppressed in comparison 
with our results, while the production of bosons could be quite significantly enhanced. The effect 
of quantum statistics is expected to give rise to a resonance amplification of the boson production 
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and will be considered elsewhere. Finally we assume, that A4 in the potential is so small that we 
can ignore self-production of (^-quanta by the field ip c . 

The case of a strong external field is considerably more complicated. There is no closed expression 
for the Green's function of the produced particles in a strong external field and no closed equation 
for (p c (t) can be written down. To overcome this problem one may use an adiabatic approach. 
In the case of the bilinear coupling to the produced particles ([53]) the calculation of the particle 
production rate can be reduced to the solution of an ordinary differential equation with time- 
dependent frequency (see e.g. GJ). It is all done in the "fixed" field approximation i.e. without 
field decay due to particle production. The corresponding decrease of the field amplitude can be 
taken into account adiabatically by imposing conservation of the total energy. In the historically 
first papers O, O, 14 1 on particle production by the inflaton field, the process was considered 



perturbatively in a spirit close to the present work, but without inclusion of the effect of back- 
reaction of the produced particles on the evolution of the inflaton field. 

The non-perturbative calculations, valid for an arbitrary strongly coupled field, p c (t), were first 
undertaken in refs. |l^, [l(|. In both papers the possibility of an enhanced production of bosons due 
to parametric resonance excitation was noticed, but it was argued that parametric resonance 
is not operative, in the framework of the model considered there, because of the cosmological 
expansion and rescattering of the produced particles. It has later been found that the resonance 
may be wide and the effects of particles leaving the resonance mode is not so important, permitting 
parametric enhancement of boson production by the inflaton decay (see ref. |17j|). 



Non-perturbative production of fermions was considered in ref. |15| for the usual Yukawa coupling to 
fermions, gtpvpip, and for an arbitrary time dependence of the scalar field, <p c (t). Concrete examples 
were given for the case of harmonic oscillations of <p, which take place if (p "live" in the potential, 
V(tp) = m 2 </? 2 /2, and the back reaction of particle production on (^-evolution is neglected. It was 
shown there that in the weak field limit the rate of particle production is described by a constant 

decay width, T, in accordance with perturbation theory. In the limit of a strong field the production 

1/2 

probability goes down as 1/V , where (fo is the amplitude of the oscillations of the field ip(t). This 
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suppression is related to an increase of the effective mass of the produced fermions, m e ff ~ g<p c , 
while the frequency of the oscillations of <p remains constant. This effect, of course, cannot be 
traced in the present work, where the weak field approximation is used. 



In a recent paper [18] a non-peturbative approach has been applied to fermion production by the 
field ip c {t) with the potential m 2 ip 2 /2 + \ip /4!. A system of coupled equations for the evolution of 
the field, (p c , and the mode functions of the produced fermions was solved numerically in one loop 
approximation. This permitted to take the back reaction of particle production on (/^-evolution 



into account. However, our equation differs from that derived in ref. [18] in the weak field limit 
(linear in the field (p c ). This difference can possibly be ascribed to a difference in our approaches 
to renormalization and the initial time singularity. As is seen, from the arguments presented here, 
we treat the initial time singularity as a usual ultraviolet singularity which can be regularized by 



the standard renormalization procedure. In ref. [18] this singularity was avoided with the help of 



a Bogoliubov transformation. In the limit of a large field strength, such that the effective mass of 
fermions, gip c , is larger than the mass of the field ip c , the calculations of ref. |18[| showed a strong 



suppression of particle production in agreement with ref. [15]. However, it is found there that the 



suppression persists even in the case of gip c <C m^,, though not for extremely small gip c . This result 



is in contradiction with the calculations of ref. [15], according to which the production rate in this 



region of parameters is well described by perturbation theory. 



Recently the approach of the paper 15] was repeated in ref. [19] for the particular case of particle 
production by the scalar field in the potential A(/> 4 . It was found that the suppression of the particle 
production at large values of the field amplitude for the case of harmonic potential, as observed 
in ]T5"|, does not take place for V = A</> 4 /4, and the production of fermions may be quite efficient. 



This conclusion was also made (a little earlier) in ref. [18]. Possible explanations of the deviation 
of these results, from the results of ref. flgl , is an influence of the A(/j 4 -term in the potential and 
corresponding excitation of higher modes, as well as an increase of the oscillation frequency of ip c 
due to a bigger amplitude of ip. 

Returning to the results of the present paper, we repeat that they are valid in the weak field 
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approximation and, as such, can be applied to the universe heating by the inflaton decay at the 
final stage of the inflationary process, when the amplitude of the inflaton field is sufficiently small. 
Even though the inflaton energy during this period only contributes a relatively small fraction of 
the initial inflaton energy, its role in the creation of the matter in the universe could be much 
more significant. The reason for that is that at the final stage of the decay, when the potential is 
approximately m 2 ip 2 /2, the inflaton behaves as non-relativistic matter. Therefore its energy density 
is amplified, with respect to the energy of the relativistic particles produced earlier, by the red- 
shift factor a(t)/ai n (t). This makes the role of the final part of the inflaton decay correspondingly 
enhanced and the temperature at the end of the decay becomes close to the lowest order perturbative 



result |20[ |. Another mechanism of suppression of particle production was discussed in ref. [21], 
where it was shown that parametric amplification might be effectively suppressed by the final state 
interaction of the produced particles. 

Another possible application of our results is spontaneous baryogenesis. The evolution of the spa- 
tially homogeneous pseudo-goldstone field, 9, which is related to spontaneous breaking of baryonic 
charge conservation, is described by the equation: 

f 2 + U'{9) = j t J% , (54) 

where J% is the baryonic charge density and / is the scale of spontaneous symmetry breaking. 
Since, as we have shown, the impact of particle production on the evolution of <p c is not given by 
the term Tip, it means that the naive identification f 2 T9 = dJ%/dt is not correct. Correspondingly 
the baryon asymmetry is not given by 5Jg ~ Y59 (for a more detailed discussion see refs. [||, ||). 
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Figure Captions: 

Fig. 1 Trilinear scalar coupling, ftpx 2 -, and harmonic potential, V = m 2 (p 2 /2. The decrease of 
the amplitude is given by the envelope function, exp(— aTt/2), T = f 2 /8-irm. Fig. la: m = 1,/ = 
0.7, A = 0, a = 1.03. Fig. lb: m = 1, / = 1.4, A = 0, a = 1.25. 

Fig. 2 Trilinear scalar coupling, f(fx 2 , and potential, V = X(f) A /4. The decrease of the amplitude 
is given by the envelope function, 1 — 0.016a/ 2 (t — ti n ). Fig. 2a: m r = 0, / = 0.5, A = 1, a = 0.94. 
Fig. 2b: m r = 0, / = 1.0, A = 1, a = 0.92. 

Fig. 3 Quartic scalar coupling, \2if 2 x 2 , and harmonic potential, V = m 2 (j) 2 /2. The decrease 
of the amplitude is given by the envelope function, (po (1 + aX 2 ,/ (16tt)ipq /m(t — ti n )) 1 ^ 2 . Fig. 3a: 
m = 1, A2 = 0.5, A4 = 0, a = 1. Fig. 3b: m = 1, X 2 = 1.5, A4 = 0, a = 1.5. 

Fig. 4 Quartic scalar coupling, A2^ 2 X 2 > and potential, 1/ = A4^ 4 /4. The decrease of the 
amplitude is given by the envelope function, y>o/(l + O.OSbaX^/ X4(po(t — U n )). Fig. 4a: m r = 
0, A 2 = 0.5, A 4 = 1, a. = 1.05. Fig. 4b: m r = 0, A 2 = 1.0, A 4 = 1, a = 1.3. 

Fig. 5 Trilinear fermion coupling, gipiptp, and harmonic potential, V = m 2 <j) 2 /2. The decrease 
of the amplitude is given by the envelope function, exp(— aTt/2), T = g 2 m/(8ir). Fig. 5a: m = 
1,0 = 0.5, A = 0, a = 1.0. Fig. 5b: m=l,g = 1.0, A = 0, a = 0.8. 

Fig. 6 Trilinear fermion coupling, g^tpf: and potential, V = X^* /A. The decrease of the 
amplitude is given by the envelope function, (fo/(l + 0.01 lag 2 -v/A^o^ — tin))- Fig. 6a: m = 0, g = 
0.5, A = 1, a = 1.0. Fig. 6b: m = 0, g = 1.0, A = 1, a = 0.82. 

Fig. 7 Quartic fermion coupling, K(p 2/ if;ip, and harmonic potential, F = m?(j) 2 /2. The decrease of 
the amplitude is given by the envelope function, [l + aipQK 2 mt/(47r)] . Fig. 7a: m = l,n = 
0.5, A = 0, a = 1.0. Fig. 7b: m = 1, k = 1.0, A = 0, a = 1.0. 

Fig. 8 Quartic fermion coupling, Kip 2 -ipif), and potential, V = X^cffi/A. The decrease of the 
amplitude is given by the envelope function, (1 + 0.075a\/A4K 2 (/?gt) 1 ^ 3 . Fig. 8a: m = 0, k = 
0.5, A = 1, a = 1.0. Fig. 8b: m = 0, k = 1.0, A = 1, a = 0.9. 
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